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ABSTRACT 

We present the analysis of the luminosity function of a large sample of galaxy clusters 
from the Northern Sky Optical Cluster Survey, using latest data from the Sloan Digital 
Sky Survey. Our global luminosity function (down to Af,. ^ —16) does not show the 
presence of an "upturn" at faint magnitudes, while we do observe a strong dependence 
of its shape on both richness and cluster-centric radius, with a brightening of M* and 
an increase of the dwarf to giant ratio with richness, indicating that more massive 
systems are more efficient in creating/retaining a population of dwarf satellites. This 
is observed both within physical (0.5 R200) and fixed (0.5 Mpc) apertures, suggesting 
that the trend is either due to a global effect, operating at all scales, or to a local one 
but operating on even smaller scales. We further observe a decrease of the relative 
number of dwarf galaxies towards the cluster center; this is most probably due to tidal 
collisions or collisional disruption of the dwarfs since merging processes are inhibited 
by the high velocity dispersions in cluster cores and, furthermore, we do not observe 
a strong dependence of the bright end on the environment. 

We find indication that the dwarf to giant ratio decreases with increasing redshift, 
within 0.07 ^ z < 0.2. We also measure a trend for stronger suppression of faint 
galaxies (below M* +2) with increasing redshift in poor systems, with respect to more 
massive ones, indicating that the evolutionary stage of less massive galaxies depends 
more critically on the environment. 

Finally we point out that the luminosity function is far from universal; hence the 
uncertainties introduced by the different methods used to build a composite function 
may partially explain the variety of faint-end slopes reported in the literature as well 
as, in some cases, the presence of a faint-end upturn. 

Key words: Large Scale Structure of Universe - Galaxies: clusters: general - Galax- 
ies: evolution - Galaxies: luminosity function, mass function - Galaxies: statistics 



1 INTRODUCTION 

Due to its integrated nature, the luminosity function (LF) 
of a specific class of objects can be used to study the distri- 
bution of luminous matter in the Universe, after taking into 
account the systematic uncertainties due to cosmic variance 
(e.g. Binggeli et al. 1988; Blanton et al. 2001; Robertson 
2010). In particular, the main focus of recent work has been 
the connection between galaxies and dark matter halos - 
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constraining various physical mechanisms governing the for- 
mation and evolution of galaxies (e.g. gas cooling, star for- 
mation, etc.). The study of the halo occupation distribution 
linked to the LF became a key factor in not only under- 
standing physical processes shaping galaxies (e.g. Peacock & 
Smith 2000; Berhnd & Weinberg 2002; Bullock et al. 2002; 
Scranton 2002) but also in providing constraints on cosmo- 
logical models (e.g. Zheng & Weinberg 2007). 

Despite the apparent simplicity of deriving the LF of 
galaxy clusters, and the many works published in the last 
few years (Lin et al. 1996; De Propris et al. 2003; Andreon 
et al. 2005; Popesso et al. 2005; Hansen et al. 2005; Gonzalez 
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et al. 2005; Zandivarez et al. 2006; Hansen et al. 2009), our 
ability to properly characterize the luminosity distribution 
of these systems has been hampered by the need to estab- 
lish well defined, statistically large and robust samples, and 
to properly combine them to address an intrinsically rimlti- 
parametric problem. In fact, the dependence of the slope of 
faint-end of the LF (i.e. the giant-to-dwarf galaxy ratio) on 
environmental and evolutionary parameters is still debated, 
as is the presence of an upturn at faint magnitudes (Hansen 
et al. 2005; Zucca et al. 2009). 

In the past, LF studies were based on selection in a 
single waveband, but today this has changed dramatically 
with surveys spanning from the ultraviolet to the infrared 
and radio bands. In particular, in the optical, large photo- 
metric surveys are now available allowing the use of cluster 
richness (a proxy for mass, see for instance Hilbert & White 
2010; Mandelbaum et al. 2010) to characterize galaxy sys- 
tems. Richness is often measured as the number of galaxies 
within a given luminosity range and within a certain dis- 
tance from the cluster center (e.g. Dalton et al. 1992; Post- 
man et al. 1996; Gal et al. 2003) allowing to stack system in 
richness bins and measure the LF over a wide range of host 
halo masses (see Gladders & Yee (2005) for different rich- 
ness definition). The purpose of the present work is to ap- 
proach the problem analyzing the uncertainties introduced 
by different reduction and analysis techniques adopted in 
the literature using a large sample of galaxy clusters with 
well defined photometry. Thus we can establish a firm basis 
for determining which results are robust and which depend 
on the specific choices for cluster detection made by differ- 
ent authors (e.g. Olsen et al. 1999; Postman et al. 2002; Gal 
et al. 2009). 

This paper is organized as follows: Section 2 describes the 
cluster catalog we use, discussing the cluster properties as 
well. In Section 3 the statistical background subtraction is 
discussed in detail as this is one of the main components 
required to estimate individual cluster LPs, which are pre- 
sented in Section 4. Section 5 presents the methods used in 
this paper for composing the individual LFs, either through 
a non-parametric or parametric approach. In Section 6, we 
discuss the main results obtained here, including the depen- 
dence of the LF on environment and its redshift evolution. 
Conclusions are drawn in Section 7. 

Throughout this paper we use a cosmology with Ho = 
70 km s"^ Mpc"^ (0^ = 0.3, = 0.7). 

2 DATA 

2.1 Cluster Catalog 

The Northern Sky Optical Cluster Survey (hereafter 
NoSOCS, Gal et al. 2009) is a new, objectively defined 
catalog of galaxy clusters drawn from the Digitized Sec- 
ond Palomar Observatory Sky Survey (DPOSS). Clusters 
are detected in two steps. First, the positions of galaxies 
from DPOSS are used to generate adaptive kernel density 
maps. Then, S-Extractor is run to detect peaks in the den- 
sity maps, which are identified as cluster candidates. A de- 
tailed description of the survey and cluster detection tech- 
nique can be found in Gal et al. (2000, 2003). Details of the 
photometric calibration and star/galaxy separation are dis- 
cussed in Gal et al. (2004) and Odewahn et al. (2004). The 



original catalog has been recently updated by improving the 
definition of bad areas, masking out very bright objects on 
the original DPOSS data, and by performing photometric 
redshift and cluster richness (Ngai) estimates for all detected 
clusters (Gal et al. 2009). Enhancement in the photometric 
redshift measurements based on DPOSS photometry involve 
measurement of the fore- and background galaxy contaminar 
tion in the cluster area. The 3-cIipped medians of the color 
and magnitude distributions from ten background regions 
arc used as the background correction for that cluster. The 
redshift estimator is run ten times for each cluster candidate. 
During the photometric redshift measurement process, the 
cluster positions arc also recomputed, leading to a ^ 30% 
improvement over the Zphot ~ 0.033 of Gal et al. (2003). At 
each iteration of the zphot computation, the median posi- 
tion of the galaxies within a 1 h"'^ Mpc radius of the pre- 
viously determined center is calculated, and this is taken as 
the new cluster ccntroid for the next iteration of the photo- 
metric redshift estimation. The resulting final NoSOCS cat- 
alog consists of 15, 502 clusters at redshift z< 0.4. The cat- 
alog comprises two separated areas of the sky: i) the North 
Galactic Pole (NGP) region covering 8494 deg^, and ii) the 
Southern Galactic Pole (SGP), corresponding to a 2917 dcg'^ 
sky region. The total contamination of the NoSOCS cluster 
sample amounts to about 8% (for details see Gal et al. 2009). 
For very rich clusters {Ngai > 50), contamination is negli- 
gible, while it rises above 5% for Ngai < 20 — 25. In the 
present work, to keep contamination rate below 5%, we fo- 
cus on a sub-sample of NoSOCS, selecting only clusters with 
Ngai > 25, in the redshift range 0.07 ^ z < 0.2 (as outside 
these limits the sample is poorly defined; see Gal et al. 2009). 

Since the LFs are computed using galaxy photometry 
from the Sloan Digital Sky Survey (SDSS), hence taking 
advantage of its better photometric accuracy relative to 
DPOSS (see below), our analysis is restricted to clusters 
found in the area imaged by SDSS. This is done by requir- 
ing that, for a given cluster, the entire region we use to 
derive the LF (see § 3) is enclosed inside the SDSS area. 
These selections result into a final sample of 1, 451 galaxy 
groups and clusters, whose distribution on the sky is plotted 
in Fig. 1 

2.2 Galaxy Catalog 

The galaxy catalog is obtained from the Sloan Digital Sky 
Survey^ 6th data release (hereafter SDSS-DR6), covering a 
total sky area of 9583 deg^ (Fukugita et al. 1996; Gunn et al. 
1998; York et al. 2000; Adelman-McCarthy et al. 2008). For 
each cluster, we select all galaxies within a 10 x 10 Mpc'^ 
region, centered on the cluster centroid. The SDSS photom- 
etry of point-like sources is 95% complete down to an r'-band 
model magnitude of mr = 22 (Stoughton et al. 2002). Since 
the SDSS star/galaxy classification is still reliable down to 
nir ^ 21.5 (Lupton et al. 2001: see also Capozzi et al. 2009), 
we select all galaxies down to the latter magnitude limit, 
and adopt this value as the apparent completeness magni- 
tude m-c of the galaxy catalog of a given cluster. This choice 
also ensures that there are no threshold effect when apply- 
ing the K-corrections discussed below, since the complete- 

^ http://www.sdss.org/ 
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Figure 1. Projected distribution of the sub-sample of NoSOCS 

clusters analyzed in the present work (gray dots). Black solid 
circles show the random regions where the global background is 
computed, to check robustness of the local background determi- 
nation (see § 3). 



ness limit of the galaxy catalog is 0.5 mag deeper than our 
adopted cut. The completeness absolute magnitude of the 
NoSOCS sub-sample used in this work, then ranges from 
around Mc = — 16 for the low-redshift clusters {z ~ 0.07), to 
about Mc = —19 for the upper redshift limit of x ~ 0.2. We 
retrieve only galaxies with clean photometry from SDSS, by 
selecting only PRIMARY objects with SDSS photometry 
flags ^ set following Yasuda et al. (2001). In order to consider 
only regions for which the galaxy catalog has homogeneous 
photometric accuracy and completeness characteristics, we 
also exclude objects within circular regions around bright 
stars and large galaxies, from the Tycho-2 and RC3 cata- 
logs (de Vaucouleurs et al. 1991), respectively. The radius r 
of the masked regions is chosen following the prescriptions 
of Gal et al. (2009). For Tycho-2 stars we set r = 2' for 
mxycho < 7.0, r = 1.5' for 7.0 < mxycho ^ 8.0, r = 1.0' for 
8.0 ^ mTycho ^ 9.5, while for galaxies in the RC3 catalog 
we adopt r = 5 x racs for racs < 25" and r = 8 x rac3 for 
rac3 > 25". 

The galaxy LF is measured in the r' band, using model 
galaxy magnitudes from the SDSS Photo pipeline (Lupton 
et al. 2001; Stoughton et al. 2002). Magnitudes arc corrected 
for Galactic extinction according to Schlegel et al. (1998). 
For all galaxies in the region of a given cluster, apparent 
magnitudes are converted to absolute magnitudes by the 
relation: 

M = m-51ogio(I>L/10pc)-K(z), (1) 

where Dl(z) is the luminosity distance, and K(2:) is the k- 
correction. In order to compute Dl(2), we assume all galax- 
ies in the given region to be at the same redshift as the 
cluster. While this is correct for cluster galaxies (considering 
the lower redshift limit of the NoSOCS catalog), it is cer- 
tainly incorrect for foreground/background galaxies. How- 

^ See also 

http : //cas . sdss . org/ astrodr6/en/help/docs/realquery . asp#e 



ever, the contribution of field galaxies is statically removed 
when computing the LF, making the computation of T)i^{z), 
on average, statistically correct. 

In order to test the robustness of our results with re- 
spect to the method used to calculate the K(z) , wo adopt two 
independent approaches, by (i) using an average k-correction 
for all galaxies, and (ii) estimating a specific k-correction 
for each galaxy. In case (i), we adopt the k-correction for 
elliptical galaxies from Fukugita et al. (1995) at the cluster 
redshift. In case (ii), wc estimate the V^{z) by using the soft- 
ware kcorrect (version 4_1_4, Blanton et al. 2003), through 
a rest-frame filter obtained by blue-shifting the throughput 
curve of the SDSS r-band by a factor (1 + zo). For zo — 0, 
one recovers the usual k-corrcction. As in previous works 
based on SDSS data (e.g. Hogg et al. 2004), we have adopted 
zo = 0.1. For galaxies at redshift z = zo, the k-correction 
is equal to — 2.51og(l -|- zo), independent of the filter and 
the galaxy spectral type. Hence, since the value of zo = 0.1 
is very close to the median redshift of the NoSOCS clus- 
ter catalog (see Gal et al. 2004), this choice of zo allows 
uncertainties on k-corrections to be minimized (see Blanton 
et al. 2003). We run kcorrect using the ugriz SDSS model 
magnitudes, and the best estimate available for the redshift 
z of each galaxy, i.e. either the spectroscopic redshift or the 
photometric estimate, when the former is not available. The 
two approaches have different advantages and drawbacks. In 
case (i) , we are assuming that early- type spectral types dom- 
inate the cluster galaxy population at all magnitudes. While 
this is only a rough approximation, it allows us to avoid 
bringing the uncertainties on the k-correction of each single 
galaxy into the computation of the LF. On the other hand, 
method (ii) implies a larger uncertainty on the K{zys, but 
corrects each galaxy according to its proper spectral type, 
inferred from the photometric information. Since the fore- 
ground/background contaminants are statistically removed 
from the LF, both methods shoxild be statistically insen- 
sitive to the K{z) of field galaxies. Fig. 2 shows the LFs 
of galaxies obtained with both methods to estimate the k- 
correction, for two richness bins of the parent clusters. The 
two methods provide very similar LFs in both cases. In what 
follows, all results are obtained by applying method (i), but 
all of our results remain essentially unchanged when using 
method (ii). 

2.3 Cluster Properties 

In order to analyze the environmental dependence of the 
LF, we derive it as a function of the richness of the parent 
clusters and the cluster-centric distance, obtained as follows. 
The optical richness of a cluster of galaxies, i.e. the number 
of galaxies in a given magnitude range within a given phys- 
ical region of the cluster, is a proxy for its mass (Kravtsov 
et al. 2004; Gonzalez et al. 2005; Popesso et al. 2007; Hilbert 
& White 2010; Mandelbaum et al. 2010). As there is no best, 
objective prescription to measure the richness parameter, we 
have performed the analysis by using two different richness 
estimates. 

i) For each NoSOCS cluster. Gal et al. (2009) measured 
the richness parameter, Ngai-, as the background subtracted 
number of galaxies in the cluster, within an aperture of 
[§.5 Mpc, in the (r-band) magnitude range of —22 to —19. To 
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Figure 2. Best-fits LFs obtained using the the Maximum- 
Likehhood approach (Sec. 5.2) by applying an average k- 
correction value from Pukugita et al. (1995) (solid line) or k- 
corrccting each single galaxy independently with the software 
kcorrect (dashed line). Gray and black colors show the cases of 
rich and poor clusters, classified according to the richness param- 
eter Richn.gDgg (see Sec. 2.3). 



take advantage of the more accurate SDSS photometry (rel- 
ative to DPOSS), we re- measured cluster richness according 
to this same definition using SDSS photometry. Hereafter, 
this updated richness parameter is indicated as Richn. sdss^- 
This approach has the advantage of being simple and com- 
monly used in the literature, as well as providing richness 
estimates well correlated to cluster mass (Lopes et al. 2006). 
However, it becomes more and more uncertain for poorer 
systems, as it uses only a limited range of the magnitude 
distribution of cluster galaxies. Notice also that the lumi- 
nosity range in the definition above extends two magnitudes 
below M* , possibly introducing a spurious dependence of 
the LF faint-end shape on cluster richness. 

ii) To minimize the drawbacks in the definition of 
Richn. SDSS, we also define a new richness estimate, 
Richn. ML, as the integral of the best-fitting Schcchtcr func- 
tion to the cluster LF in the luminosity range Mr ^ —21.0. In 
order to account for the possible dependence of the Schechter 
fit on cluster richness (see § 6.1 ), we estimate Richn.ML us- 
ing an iterative procedure. First, we split our cluster sample 
according to Richn. sdss, and derive the Schcchtcr fit to the 
LF for each richness bin. For each cluster in the bin, the 
LF fit is re-scaled to match the number counts of galaxies 
brighter than —21 . This provides a first richness estimate, 
Richn.ML- We then re-arrange the cluster sample based on 
Richn.ML, Btnd repeat the procedure, obtaining the second, 
final richness estimate, Richn.ML- The Schechter function 
fits are obtained with the Maximum-Likelihood approach 
described in Sec. 5.2. Notice that this iterative procedure is 
preferable to measuring the Richn.ML from the Schechter fit 
to the LF of single clusters, since results of single fits exhibit 
a large measurement scatter (see §6). Richn.ML is designed 



^ We note that Gal et al. (2009) actually used an iterative method 
in order to include k-corrections in the richness estimates. This 
refinement is not implemented in our procedure. 
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Figure 3. Comparison of Richn.gnss and Richn.ML- Notice the 
good correlation among the two estimates. 



to probe mainly the LF normalization at M* , i.e. the abun- 
dance of giant galaxies in the cluster, independent of the LF 
faint-end slope. On the other hand, it is based on the in- 
tegral of a parametric fit re-scaled to the whole magnitude 
distribution of galaxies in a cluster, hence being virtually less 
affected by Poissonian noise on number counts with respect 
to richness estimates obtained by the number of galaxies in 
a given magnitude range. To further verify the dependence 
of Richn.ML on the LF faint-end slope, we simulated a set 
of clusters with Richn.ML ranging from 4 up to 15, i.e. the 
range covered by our cluster sample, and intrinsic a = —0.8. 
We then recomputed Richn.ML, fitting the LF with a forced 
to be larger/smaller by 0.6 than its best-fit value, i.e. > la in 
the large majority of cases of Table 2; this would correspond 
to a situation where a cluster is assigned to a completely er- 
roneous group and thus its Richn.ML is measured using the 
wrong Schechter model. Even is such extreme scenario the 
variations in Richn.ML arc less than 16% in all cases, and 
anyway always below the poissonian uncertainties affecting 
richness estimates based only on galaxy counts. 

Fig. 3 compares Richn.sDSS with Richn.ML- Despite the 
different definitions, a good correlation is observed between 
the two sets of measurements. More generally, we verified 
that our results remain unchanged regardless of which rich- 
ness estimate is used. In particular, the dependence of the 
LF on environment is the same using either Richn.ML or 
Richn. SDSS, even though in the latter case the trends re- 
ported in §6.1 appear somewhat weaker due to the rea- 
sons discussed above. For brevity, throughout the rest of 
the work, we will only show results obtained for Richn.ML- 

In order to estimate the cluster-centric distance of 
galaxies in different clusters, wo use both fixed and, when 
available, characteristic radii -R200 - The characteristic radius 
of a cluster, R200, is the radius within which the mean in- 
ner density is 200 times the critical density, Pc{z), of the 
Universe at the cluster redshift. N-body simulations suggest 
that the bulk of the virializod mass of a cluster is gener- 
ally contained within this radius (e.g. Carlberg et al. 1997). 
Values of R200 were computed for a sub-sample of NoSOCS 
clusters from Gal et al. (2009), assuming that the radial 
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distribution of galaxies within a cluster follows that of the 
dark matter, and neglecting possible variations of the mean 
mass of galaxies with environment. Of the 1, 451 clusters an- 
alyzed in the present work, R200 values are available for 814 
clusters. For this sub-sample, the LFs are derived in 

• three circular regions, with outer radii of 0.2 R200, 
0.5 R200 and R200, all centered on the cluster centroids; 

• an outer annulus, with inner and outer radii of 0.5 and 

1 Ji200. 

Deriving the LFs within regions sampling different frac- 
tions of a dynamical radius, such as Ji2oo, rather than 
fixed-size apertures, can actually provide a more physically- 
meaningful way to compare galaxy populations at different 
cluster-centric radii. However, in order to fully exploit the 
extensive statistical power of the NoSOCS cluster sample, 
we also perform a fixed-aperture analysis of the LF using the 
entire sample of 1,451 clusters. This also allows us to per- 
form a more direct comparison with previous works, where 
R200 measurements were not available. For the entire sam- 
ple, wc derive the LFs within five fixed-size apertures: 

• three circular regions of radii: 0.5, 1.5 and 3.0 Mpc; 

• two concentric annuli, with 0.5 < r ^ 1.5 Mpc, and 
1.5 < r < 3.0 Mpc. 

Throughout this paper, results using fixed apertures al- 
ways refer to the entire cluster sample, while results within 
the physical apertures are obtained using the sub-sample of 
814 clusters. 



3 BACKGROUND STATISTICAL 
SUBTRACTION 

The LF of a galaxy cluster is defined as the number of galax- 
ies in the cluster as a function of luminosity. The primary dif- 
ficulty in measuring the LF is that of assessing the member- 
ship of galaxies along the line of sight. Furthermore, projec- 
tion effects tend to mimic the presence of a large population 
of dwarf galaxies, hence producing steep faint slopes (Val- 
otto et al. 2001). Ideally, one would need spectroscopic red- 
shifts for each individual galaxy in the cluster area, but this 
is unfeasible for a plethora of reasons. Spectroscopic mea^- 
surements axe extremely time demanding if possible at all: 
i.e. when the sample of clusters is large, when the faint end of 
the galaxy population has to be analyzed, when dealing with 
high redshift clusters. Several authors have attempted to use 
photometric redshifts to assess cluster membership but, even 
though these are effective at reducing the back/foreground 
contamination, the latter remains non-negligible and statis- 
tical corrections are still required to remove the contribution 
of contaminant sources (Tanaka et al. 2005, 2007; Rudnick 
et al. 2009; Capozzi et al. 2009). 

The statistical background subtraction is performed by 
estimating the contribution of non-cluster members to the 
number counts of galaxies in the cluster direction, by mea- 
suring the projected number counts of field galaxies outside 
the cluster region. Two approaches can be used, measuring: 
(i) the "global" density of field galaxies on large angular ar- 
eas (Gladders & Yee 2005; Hansen et al. 2005), or (ii) the 
"local" background either in control fields close to the clus- 
ter, or in annuli centered on the cluster centroid (Paolillo 
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Figure 4. Top panel: Average background counts per unit area 
as a function of r-band magnitude (filled circles); error bars 
are smaller than symbols. Results from Popesso et al. (2004) 
(empty triangles) are also shown. The solid line plots the count- 
magnitude relation predicted from Yasuda et al. (2001) in the 
magnitude range 12 ^ magr =S 17: the dotted line represents its 
extrapolation to fainter and brighter magnitudes. Bottom panels: 
percentage difference between local background counts per unit 
area, outside 3 Mpc, and the: i.) local background counts out- 
side 4 Mpc: ii.) local background counts (outside 3 Mpc) for the 
200 poorest clusters in the NoSOCS sample; iii.) local background 
counts (outside 3 Mpc) for the 200 riciicst clusters in the NoSOCS 
sample, iv.) global background (see the text). 



et al. 2001; Goto et al. 2002; Popesso et al. 2005). The dif- 
ferences between these two methods have been extensively 
analyzed in the literature, with most authors finding no sig- 
nificant difference between them (e.g. Driver et al. 1998; 
Goto et al. 2002; Hansen et al. 2005; Popesso et al. 2005; 
Barkhouse et al. 2007). On the other hand, as noted by Pao- 
lillo et al. (2001, also see below), when one wants to isolate 
just the main cluster signal without analyzing the structure 
correlated with the cluster, the local background approach 
is preferable, as it takes into account possible background 
variations in the cluster region, caused by the large-scale 
structure within which the clusters are embedded. 

Wc therefore derive the galaxy LF using a local back- 
ground approach. The local background is estimated within 
a 10 X 10 Mpc region, centered on the cluster centroid, 
outside a radius of 3 Mpc (about two times the Abell ra^ 
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dius), where the contamination from cluster galaxies is ex- 
pected to be negligible. In order to avoid over-estimating 
the background level because of the possible presence of 
back/foreground galaxy groups, we follow the approach of 
Paolillo ot al. (2001). Wc generate a density map of galaxies 
in the background region by convolving the projected dis- 
tribution of galaxies with a Gaussian kernel of cr = 250 kpc 
in the cluster rest frame (the typical size of a cluster core). 
Then, we mask out all density peaks, above the 3a level, 
from the background region. Masked out regions cover, on 
average, 2 — 3% of the whole background area and contain 
less than 2% of all background galaxies. Clusters for which 
the area of the masked regions is larger than 10% of the total 
background one are excluded from our analysis. The rmm- 
ber counts of galaxies in the remaining region, which we call 
the "local background" , is adopted to estimate the expected 
background counts in the cluster direction, i.e. one of the 
regions where the LF is derived (see Sec. 2.3). Fig. 4 (top 
panel) shows the local background counts, per unit area, av- 
eraged among all the clusters in the NoSOCS sample. Good 
agreement is found with Popesso et al. (2004), who esti- 
mated background counts within randomly selected fields, 
and with the count-magnitude relation expected for a ho- 
mogeneous galaxy distribution in a universe with Euclidean 
geometry, as obtained by Yasuda et al. (2001). 

Fig. 4 also shows some tests we performed to check the 
robustness of the local background determination. The bot- 
tom panels show the fractional difference between averaged 
local background counts obtained outside 3 Mpc and those 
measured: i) outside 4 Mpc; ii) outside 3 Mpc but only for 
the 200 poorest clusters in the NoSOCS sample; iii) outside 
3 Mpc but only for the 200 richest clusters. Panels i.-iii. 
show no appreciable difference (< 1%, on average) in num- 
ber counts in the two background regions, implying that, on 
average, we are not overestimating the background counts, 
as might be the case if some residual signal from the clus- 
ter would be still detectable in our "local" background. If 
present, such effect would indeed be less important at larger 
cluster-centric distances (panel i.), and/or produce a ficti- 
tious higher background level for richer clusters (panel iii.), 
while no measurable difference is instead observed. Sheldon 
et al. (2009), in their weak lensing analysis, measure the ex- 
cess number density due to the cluster in the background 
region to vaxy between ~ 1% for poorest groups to ~ 3.5% 
for rich clusters. Poorest groups are excluded from our clus- 
ter sample, and hence the difference in the excess number 
densities caused by poor and rich clusters in the surrounding 
background region should be lower than ~ 2.5%, measured 
over the whole magnitude range, in rough agreement with 
what found in this work. 

As a further test, we also extracted a "global" back- 
ground from 42 control fields, selected within the SDSS area, 
containing no known galaxy clusters. The control fields have 
a radius of 60' each, covering a total area of ~ 130 deg^. 
Their distribution on the sky is shown in Fig. 1 (black solid 
circles). The difference between the counts measured in this 
"global" background, and those obtained in the "local" one 
(outside 3 Mpc), is plotted in panel iv. of Fig. 4. While for 
bright magnitudes, m,. < 14, the two background estimates 
are consistent within the errors, the global background tends 
to be systematically lower than the local one for fainter mag- 
nitudes, in agreement with the findings of Paolillo et al. 



(2001). Since wc can reasonably exclude that our back- 
ground local estimate is over-estimated because of contam- 
ination from cluster galaxies (see panel i. in Fig. 4), we can 
conclude that the local background also includes the con- 
tribution of the large-scale structure around clusters and 
groups of galaxies, which is instead not accounted for by 
the number counts in the control fields. Throughout the 
present work we then compute the LFs by using the lo- 
cal rather global background determination. We note that, 
while many authors have claimed no significant diS'erenccs 
in the LF when using either the local or global background. 
Fig. 4 of Popesso et al. (2004) reveals, in agreement with our 
findings, a tendency for the LF to have a shallower faint- 
end slope when using global background counts, implying a 
differential slope between the local and global background 
galaxy counts, which could affect the slope of the faint-end 
slope of the resulting LF. 



4 INDIVIDUAL LUMINOSITY FUNCTIONS 

In order to derive the LF of galaxies in individual NoSOCS 
clusters, for each cluster we estimate the number counts of 
galaxies in the local background region (see previous section) 
and subtract them from the number counts of galaxies in the 
cluster region (see § 2.3). Number counts are computed in 
half-magnitude bins. Background counts are rescaled to the 
effective area of the cluster region, accounting for excised 
areas due to bright objects in both background and cluster 
regions (§ 2.2). Errors on the individual LFs are measured 
following Paohllo et al. (2001). 

As an example of the procedure to derive the single LFs, 
Fig. 5 exhibits background and cluster number counts, as 
well as the resulting LF for two clusters among the poorest 
(NSC14787, z = 0.132, top panels) and richest (NSC09718, 
z = 0.137, bottom panels) structures in the NoSOCS sam- 
ple. While for rich structures the LF is significantly detected 
above the background level, poor structures are affected by 
their very low density contrast, as shown by the large un- 
certainties on background-subtracted number counts at both 
the bright and faint ends of the LF. 

In order to characterize the properties of the individ- 
ual LFs, we fit them by a parametric model given by the 
Schechter function: 

$(M)dM = 100-4(M*-M)(. + 1) ^^qOA(M'-M)^ ^^^2) 

where M is the galaxy magnitude, $* is the normalization 
factor, M* is the characteristic knee magnitude and a is 
the faint-end slope of the LF. For each cluster, we fit the 
background-subtracted counts using a minimization pro- 
cedure. To calculate the x^, the LF model is integrated over 
each magnitude bin. Uncertainties on the best- fitting param- 
eters, "1?*, M* , and a, are determined by marginalizing each 
parameter over the remaining ones. Results are presented in 
Sec. 6. 



5 COMPOSITE LUMINOSITY FUNCTIONS 

To analyze the dependence of the cluster LF on different 
properties, such as cluster richness, redshift, and cluster- 
centric distance of galaxy populations, we bin the NoSOCS 
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Figure 5. Left panels. Galaxy counts extracted within R200 
(solid circles) , and in the local background (empty circles) for the 
poor cluster NSC14787 (top) and rich cluster NSC09718 (bot- 
tom), both located at z 0.13. Right panels. Background sub- 
tracted number counts, i.e. individual LP's, of the two clusters. 
The dashed line corresponds to the completeness limit of the 
SDSS photometry (Sec. 2.2); the solid line shows the best fit mod- 
els. Poor structures, with low signal LFs, might lead to unrealistic 
results in the LF best fit (as for NSC14787, top right panel), which 
cause the large spread in the LF results, shown in Fig. 9. 



sample with respect to each quantity, and derive the com- 
posite LF of galaxies in each bin. The composite LFs are 
derived using two alternative approaches: by performing a 
weighted average of the individual cluster LFs (Sec. 5.1), and 
by performing a simultaneous Maximum-Likelihood (here- 
after ML) fit to all the individual cluster LFs in the given 
bin (Sec. 5.2). The first approach is non-parametric, i.e. we 
make no prior assumptions on the shape of the LF, while 
the ML fit assumes a given functional form. 



5.1 Non-parametric Approach: 
LF 



"cumulating" the 



The cumulative LF (hereafter CLF) of a given sample of 
galaxy clusters is the weighted mean of the individual LFs. 
Previous works have focused on the differences among alter- 
native cumulation approaches, and the following two meth- 
ods have come out as the most reliable ones: 



Colless method. CoUess (1989) derived the CLF as: 



(3) 



where Ncj is the number of galaxies in the j-th bin of the 
CLF, A'ij is the number of galaxies in the j-th bin of the i-th 

cluster LF, rrij is the number of clusters contributing to the 
j-th magnitude bin, A^io is a normalization factor and: 

= ^iVio. (4) 



For the i-th cluster, the Nio is defined as the number of 
galaxies brighter than the completeness absolute magnitude, 
Mc, of the cluster sample. For a flux- limited survey, Mc co- 
incides with the completeness magnitude of the most distant 
cluster. The error on the resulting CLF is 



m: \l Nio 



(5) 



where crjVij is the statistical uncertainty on Nij (see Sec. 4). 
For the NoSOCS sample, we compute the individual LFs 
as described in Sec. 4. The completeness magnitude of the 
sample is Mc = —19 in the r' band (see Sec. 2.2), and Nio is 
hence computed as the field-corrected number of galaxies, 
of the i-th cluster, brighter than —19. 

GMA method. Garilli et al. (1999) proposed an alternative 
cumulation method, where one weights each cluster LF ac- 
cording to the number of galaxies contained within an adap- 
tive magnitude range. In this approach, the CLF is given by: 



(6) 



where A^cj and Mj arc defined as in the Colless method, mj is 
the rnimber of clu.stor.s with completeness magnitude fainter 
than the j-th bin, and Wi is the normalization factor of each 
cluster. For the i-th cluster, we denote as Mc,i its complete- 
ness (absolute) magnitude. The Wi is given by the ratio of 
the number of galaxies brighter than Mc,i in the cluster, to 
the average number of galaxies brighter than Mc,i among 
all clusters whose completeness magnitude goes fainter than 
Mc,i. The error on the CLF is 




(7) 



where aN,, is defined as for the Colless method. 



The main diflterence between the two cumulation meth- 
ods is that Colless weights each LF by the number of galax- 
ies in a fixed magnitude range, while the GMA takes advan- 
tage of the different completeness magnitude of each cluster, 
hence fully exploiting all the available information in the 
data. In principle, both methods should provide identical 
results, at least as far as all clusters exhibit, in a statisti- 
cal sense, the same LF, i.e. the cluster LF is "universal" . 
Fig. 6 compares the CLFs obtained by the two methods for 
all NoSOCS clusters within a region of radius 0.5 Mpc. The 
methods agree very well on the bright end of the LF, while 
at faint magnitudes, we sec a trend for the Colless method 
to produce a steeper faint-end than that of the GMA CLF. 
A trend in the same direction has been already reported 
by Popesso et al. (2005) (hereafter P05), who detected a dra- 
matic drop-off in the CLF faint-end when using the GMA, 
rather than Colless, method. P05 explained this effect as a 
result of a strong correlation among the GMA weights, Wi, 
and the completeness magnitude Mc,i of the clusters. This 
would actually lead to down-weighting the LFs of clusters 
having a deeper completeness magnitude hence producing a 
sharp decline in the CLF faint end. Fig. 7 plots the GMA 
weights for the NoSOCS sample as a function of the clus- 
ter completeness magnitude, Mc,i. In contrast to P05, we 
find no strong correlation between Wi and Mc,i (cfr Fig. 5 
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Figure 6. Composite LFs (within the cluster region of radius 
0.5 Mpc) obtained applying the Colless and GMA methods plot- 
ted, respectively, as black and gray solid lines. 



o 
in 



o 







GMA 





















-18.0 



-17.0 



M 



-16.0 



-15.0 



c,i 



Figure 7. Dependence of the GMA weight, Wi, on the complete- 
ness magnitude, Mc^i, of NoSOCS clusters. Notice that no rele- 
vant correlation is detected. 



of P05), implying that the difference between Colless and 
GMA CLFs is not caused by the trend found by P05 

The difference between the Colless and CMA methods 
is instead related to the fact that the GMA computes the 
weights by also counting galaxies in the faint end of the in- 
dividual LFs. In the case where the bright part of the LF is 
identical for all clusters in the sample, all of the individual 
LFs have the same weight in the Colless method, regardless 
of the shape of the faint-end slope. For the GMA, the Wj's 
are the same only if the faint-end part of the LFs is the same 
for all clusters. If this is not the case, clusters with a steeper 
faint-end have a smaller w^^ , i.e. are down- weighted in the 
CLF (see Eq. 6). In other words, if the individual cluster 
LFs have different faint-end slopes, the GMA will give more 
weight to those clusters with a shallower LF faint-end, hence 
producing the difference between the two CLFs, as observed 
in Fig. 6. The variation of the faint-end slope among clus- 
ters can be either (i) intrinsic, i.e. the LF is not universal, 
or (ii) statistical, because of the non-Poissonian background 
fluctuation on the cluster angular scale (see Sec. 3). Poisso- 
nian uncertainties (on both cluster and field counts) do not 



We notice that P05 define the Wi as the ratio of the number 
of galaxies brighter than Af,;,i in a cluster, to the number, rather 
than the average number, of galaxies brighter than Mc,i. This 
might explain the discrepancy between the lack of trend in Fig. 7 
of this paper, and the strong trend shown in Fig. 5 of P05. 



alter the shape of the individual cluster LF, as the errors on 
different magnitudes are not correlated. This is not the case 
for the non-Poissonian contribution to the error budget. As 
shown in Sec. 6, point (i) is certainly important in driving 
the difference of the CLFs from the two methods in Fig. 6, 
as the LF faint-end slope turns out to depend significantly 
on cluster richness. 

While both the GMA and the Colless methods have 
the advantage of being non-parametric, they are also signif- 
icantly affected by statistical fluctuations in the individual 
LFs, in particular when the individual LFs have low S/N 
ratio, i.e. for poor groups and for the outskirt regions of 
clusters, as well as at the extreme faint end of the LF. At 
faint magnitudes, statistical fluctuations in the background 
signal can even lead to negative values of fleld-corrected 
number counts of the individual LFs, possibly making the 
weighted mean in Eq. 3 and Eq. 6 ill-defined. The GMA 
method is more sensitive to this issue than the Colless one, 
as the weight for each LF is computed by including also the 
faint range of the LF, where background fluctuations are 
more important. Therefore, while both cumulation methods 
can be applied straightforwardly to well controlled samples 
of rich structures, caution should be taken when analyzing 
individual LFs with a variety of S/N ratios and completeness 
limits. 



5.2 Pcirametric Approach: the Maximum 
Likehhood Technique 

An alternative approach to derive the composite LF of 
a cluster sample is to perform a simultaneous Maximum 
Likelihood flt of number counts for all clusters. While this 
method has the drawback of being parametric, in that a 
given analytic functional form of the global LF has to be as- 
sumed, it has several advantages relative to the cumulative 
approach. First, the data are not binned. Second, no correc- 
tion has to be applied for incompleteness at the faint end, 
and third, it is not affected by the issue of negative values 
of field-subtracted number counts. 

The ML approach and its advantages have been thor- 
oughly described by Andreon et al. (2005); henceforth we 
shortly describe only its implementation for the analysis of 
the NoSOCS sample. First, we have to adopt a given func- 
tional form for the global LF. As noted in Sec. 4, a single 
Schcchter function provides a reasonable tool to analyze the 
LF of individual clusters, where the uncertainties on number 
counts usually do not allow a detailed analysis of the shape 
of the LF. Now instead we consider a Schechter plus a Log- 
normal function. The latter term is used to describe the LF 
of Brightest Cluster galaxies (BCGs), which are known not 
to follow the Schechter distribution typical of faint galax- 
ies (Thompson & Gregory 1993; Biviano et al. 1995; Hansen 
et al. 2009). 

To perform the ML fits, we assume that each galaxy 
is extracted from a probability distribution consisting of the 
above model plus a second order power-law, representing the 
background component. For each cluster, we first fit the sec- 
ond order power-law to the local background. The best fit 
is rescaled to the angular area of the cluster region. Galgixy 
counts in all the cluster regions are then fitted simultane- 
ously, by keeping fixed, for each cluster, its rescaled back- 
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Table 1. Results ML fit within 0.5 R200- 



a M* 



ground power- law component ^. The fitting parameters are 
the characteristic magnitude, M* and the faint-end slope, 
a, of the Schcchtcr function, the central magnitude, M^qq, 
and width, asca, of the BCG's component, and the two 
normalization factors of the Schechter and log-normal func- 
tions. 

The normalization factors are let free to vary from clus- 
ter to cluster, while the other parameters are set to have the 
same values for all clusters. The ML fits are performed using 
the L-BFGS algorithm as for the fitting of the individual LFs 
(Sec. 4). L-BFGS is well suited for optimization problems 
with a large number of dimensions, as is the case for the CLF 
fitting, because it never explicitly forms or stores the Hes- 
sian matrix (Lu et al. 1994), still allowing upper and lower 
constraints for each variable to be measured. We also con- 
strain the fitting parameters of the log-normal component 
of the CLF to be in the range of -23.5 < Mgca =S -22.5 
and 0.3 ^ ctbcg ^ 1-0, based on the typical range spanned 
by BCGs (e.g. cf Hansen et al. 2009). Confidence limits are 
evaluated by marginalizing over all unwanted free parame- 
ters. 



6 RESULTS AND DISCUSSION 

The measurement of the LF of galaxy clusters is an intrin- 
sically multiparametric problem. Local and global galgixy 
density, density of the intra-cluster medium, age etc. are all 
known to affect the relative balance between galaxy popula- 
tions in clusters through several effects such as merging, tidal 
stripping, harassment, ram pressure stripping and strangu- 
lation, whose interplay is not yet well understood. 

We have underlined in the previous sections that the 
assumption of a universal LF shape in samples spanning a 
large range of parameters and S/N, is not merely an approx- 
imation which hides the details of the processes at work, but 
may produce biases in the final results altering the measured 
LF parameters depending on the technique used to weight 
the individual LFs. 

Nevertheless, in order to allow a comparison with pre- 
vious studies of large cluster samples, we derived a "global" 
LF of the overall sample within 0.5 R200 (see Fig.8 and Ta- 
ble 1). Differences in cumulation techniques, extraction radii 
and richness of the samples affect the shape of the LF (see 
§ 5.1 and § 6.1), making a comparison between independent 
works a difficult task. 

In general, the M* values agree very well within the 

^ Wc tested the possibility of fitting simultaneously the local 
background witli the number counts in all the cluster regions, 
but the results turned out to be indistinguishable from the case 
where the local background is fitted independently for each clus- 
ter. 




I -1.04 -1.00 -0.96 



^56 34 32 3o ^Ts 
MaQr 

Figure 8. ML fit obtained using a Schechter plus a Log-normal 
components (solid line). The individual components are shown 
as dashed lines. Data were extracted within 0.5 R200- Confidence 

levels for the free parameters of the Schechter function are also 
plotted. Gray circles and empty diamonds represent the compos- 
ite LFs obtained using the CoUess and GMA methods, respec- 
tively. 

uncertainties among different authors, while the faint-end 
slope of the LF is more debated. Once the different pass- 
bands and cosmologies are taken into account, our M* is 
indeed in fair agreement with Crawford et al. (2009) ; Garilli 
et al. (1999); Paolillo et al. (2001), with Yang et al. (2009) for 
galaxy groups in the SDSS, and with Rudnick et al. (2009) 
for red-sequence galaxies, which dominate the bright end of 
the LF. 

For the steepness of the faint end, our results are in good 

agreement with the X-ray selected clusters by Valotto et al. 
(2004) (q — 0.9 ±0.1), and in marginal agreement with Pao- 
lillo et al. (2001) (a = l.lllJJ '®). The slightly steeper faint 
end measured by Paolillo et al. (2001) is probably an effect of 
both the larger extraction radius (r ~ 2 Mpc) and the sam- 
ple of rich Abell clusters used. Trends for flatter faint end 
slopes are instead observed by Garilli et al. (1999); Crawford 
et al. (2009) (a = -O.Sllg;?^ and a = -0.84 ±0.32, respec- 
tively). While the average value of a measured by Crawford 
et al. (2009) over the whole low-redshift cluster sample is 
consistent with our findings within the uncertainties, the dif- 
ference with Garilli et al. (1999) can be attributed to their 
smaller extraction radius (< r >~ 250 kpc). 

Hansen et al. (2005, 2009), fitting a MaxBCG selected 
sample down to Mr < —19, derive a shallower faint-end. 
The latter sample however may be skewed toward poorer 
systems with respect to NoSOCS, since when differentiating 
according to cluster richness their results are in much better 
agreement with ours (see 6.1). 

On the other hand, while their M* are typically > 
0.6 mag dimmer, their cluster finding technique might fa- 
vor systems whose LF exhibits a prominent BCG compo- 
nent, resulting in a pronounced Log- normal component at 
the bright end, which is not observed in our sample ®. We 

^ Indeed, Koester et al. (2007) have shown that maxBCG catalog 
is not biased toward bright BCGs, and bright-BCG systems do 
not have satellites with systematically fainter M*. 
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speculate that the presence of a more prominent Log-normal 
component anti-correlates with the LF characteristic magni- 
tude, in the sense that systems with a stronger Log-normal 
component have dimmer satellites. 

Popesso et al. (2005) find a steeper faint end and 
brighter M* {a = -1.33 ± 0.06, M* = -22.17 ± 0.20 
within 1 Mpc h.^QQ) than we do, but their results strongly 
vary depending on the extraction region that they use. 
Their two-component fit within both 1.0 Mpc h^QQ and 

2.0 Mpc hj~QQ is in fact roughly consistent with our a 
value (a = —1.05 ± 0.13 and a = —1.03 ± 0.14, respec- 
tively), while they observe a much shallower trend within 
1.5 Mpc h^,{, (a = -0.76 ± 0.13). When comparing LFs 
extracted within the same physical radius R2()(), the steep- 
ness of the best-fit bright Schochter component from a later 
work by the same authors (Popesso et al. 2006) is in- 
stead in agreement with our findings (a = — 1.09 ± 0.09, 
while wo measure a — — 1.15 ± 0.02 over the whole R200 
sub-sample); their characteristic magnitude is only slightly 
brighter than our measured value (M* = —21.71 ± 0.16 
against M* = —21.43 ± 0.6). At the faintest magnitudes 
{Mr ^ —18) Popesso et al. (2006) detect a significant up- 
turn which is not soon in our data (see Fig. 6). However, 
we do not probe magnitudes fainter that Mr < —16; also 
note that these authors adopted the CoUess cumulation ap- 
proach which can result in a steeper LF, in particular when 
the sample is weighted using only galaxies much brighter 
than the cumulative LF limit (§ 5.1). 

Valotto et al. (2004) find that artificially steep faint 
end slopes might be caused by projection effects resulting 
from background galaxies, which cannot be corrected for by 
subtracting background fields for 2d selected clusters with 
no significant 3d couutorpart. Despite the 2d selection algo- 
rithm applied to compile the NoSOCS catalogue, the reality 
of NoSOCS clusters was verified both by photometric red- 
shifts and by X-ray analysis of a sub-sample of NoSOCS 
clusters (Lopes et al. 2006). In any case, we do not observe 
the steep LF faint end described by Valotto et al. (2004). 

In Fig. 9 we plot the results of the individual best-fits to 
all clusters within 0.5 R200. The individual LFs show a wide 
scatter, due to the combined effects of variable S/N levels 
(un-physical results arc obtained for many low S/N ratio 
or very bright completeness limit systems) and intrinsically 
different galaxy populations (see § 6.1). Despite this, the 
most likely values within the whole sample are in excellent 
agreement with the best fit of the ML fit obtained using a 
single Schechter function (plotted as a white circle). 

In the following sections we split our sample in sub- 
samples of clusters spanning small ranges of richness, 
cluster-centric distance and redshift in order to minimize cu- 
mulation biases and understand the main parameters driv- 
ing galaxy evolution in different environments. 

6.1 LF Dependence on the Environment 

Discordant results have been reported in literature in the 
past years, about M* being brighter in clusters than in the 
field, with an increasing trend for higher mass-systems (De 
Propris et al. 2003; Hansen et al. 2009; Barkhousc et al. 
2007, 2009; Crawford et al. 2009). This is in agreement 
with hierarchical models for galaxy formation and evolution, 
where the frequency of mergers increases in intermediate and 
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Figure 9. Results of the individual fit to a single Schechter func- 
tion for all clusters in our sample within 0.5 R200. Rich clusters 
(RichnML ^ 9) are plotted as gray crosses; all other clusters are 
shown as black dots. Over-plotted are the contours including 25%, 
50% and 75% of the whole distribution. A white dot shows the 
result of the ML best fit obtained using a single Schechter func- 
tion. 

high mass systems, such as rich groups and clusters, causing 
galaxies in structures to be typically brighter than in the 
field, hence resulting in a brighter M*. 

Substantial differences in the shape of the LF at dif- 
ferent selection radii have been found for clusters at both 
low (Barkhouse et al. 2007; Popesso et al. 2006; Lobo et al. 
1997; Robotham et al. 2010) and moderate-to-high red- 
shift (Crawford et al. 2009). While the slopes of the LFs are 
similar at bright magnitudes, the main differences arise at 
the faint end where the influence of the dwarf galaxies causes 
an increase in the steepness of a with cluster-centric dis- 
tance. The sampling depth and the effective cluster-centric 
distance thus have a great influence on the measured shape 
of the LF since the inclusion of different fractions of the 
dwarf galEixy population will directly impact the slope of 
the faint end. Similar results have also been reported for 
field galaxies (Xia et al. 2006). On the other hand, a lack of 
significant trends with cluster-centric radius, has also been 
reported (Hansen et al. 2009) although limited to brighter 
magnitudes (Rudnick et al. 2009). A dependence of the faint- 
end slope on the mass of the cluster has also been reported, 
suggesting that more massive clusters exhibit higher dwarf 
to giant ratios than less massive ones (De Lucia et al. 2004; 
Zandivarcz ct al. 2006; Gilbank et al. 2008). 

We start by splitting our sample into five equally pop- 
ulated sub-sets of fixed richness ranges. The top left panels 
of Figs. 10 and 11 show the ML fit for the LFs in all rich- 
ness bins, within a projected radius of 0.5Mpc and O.5R200, 
respectively. Both panels reveal a systematic change of the 
overall population across the five richness bins, in the sense 
that the lowest multiplicity bins show strong dwarf suppres- 
sion and fainter M*. We note that the decreasing relative 
normalization of the LFs is real, a result of the decreasing 
cluster richness. 

To explore the influence of the extraction region and 
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of the cluster-centric distance we fit the LPs within several 
extraction radii (both fixed and physical) for each richness 
sub-sample: an outer annulus (0.5 R200 ^ R < R200 and 
0.5 Mpc < R < 1.5 Mpc, for the physical and fixed aper- 
tures, respectively) and a larger projected radius (R ^ R2()() 
and R ^ 1.5 Mpc). Results of the ML fits are listed in Table 2 
and shown in the bottom and top right panels in Figs. 10 
and 11. 

At larger cluster-centric radii (top right panels) we ob- 
serve that the faint end slope becomes systematically steeper 
for all richness bins, respect to what observed within smaller 
apertures (0.5Mpc and O.5R200), while no dramatic change 
is observed on the bright side. Furthermore, the dependence 
on richness both of the steepness of the faint end and of 
the characteristic magnitude, is very diluted. Eventually, 
the bottom panels show that in the outer annuli clusters 
maintain a similar shape, within their 3(t errors, among all 
richness bins. This reveals that the sharp steepening of the 
faint-end and the M* brightening with richness, observed 
out to large projected radii, is mainly due to the galaxies lo- 
cated within the central cluster regions. Outside ~ 0.5 R200 
cluster galaxies appear to share similar LFs, regardless of 
the mass of the parent halo mass. It is thus likely that stud- 
ies using large apertures, either physical or fixed, are unable 
to identify such trends due to the role played by the central 
cluster regions in determining the LF shape. Similar results 
are observed at all scales: from galaxy groups (Robotham 
et al. 2010) to rich clusters (Popcsso et al. 2006). 

The decrease of the relative number of dwarf gala:xies 
towards the cluster center is most probably not due to re- 
cent merging processes, since these are inhibited by the high 
velocity dispersions in cluster cores. Also, we observe that 
the shape of the bright end docs not strongly depend on en- 
vironment; bright early- type cannot hence be the product of 
cluster environment. The most likely explanations for these 
massive galaxies is tidal collisions or coUisional disruption 
of the dwarfs, which most probably cnds-up contributing to 
the intra-cluster diffuse light. 

We point out that the normalization of the Log-normal 
component is poorly constrained in our fits especially for 
poor systems, since it is individually measured for each clus- 
ter. A more detailed analysis of the bright galaxy population 
is thus deferred to a future work. 

6.2 Redshift Evolution 

The NoSOCS sample studied in this work spans only a lim- 
ited range of redshifts (0.07 ^ ^; < 0.2). In order to study 

the evolution of the LF within this limited redshift range, we 
thus need to control the effects due to richness and extrac- 
tion radii discussed in the previous section, which may oth- 
erwise dominate possible LF variations. We hence divide our 
sample into rich (Richn.ML ^ 6) and poor (Richn.ML < 6) 
clusters in order to minimize the intrinsic variation of the LF 
shape (see § 6.1). We further limit the analysis to the cen- 
tral 0.5 Mpc, where we have observed the strongest depen- 
dence of the galaxy population with the environment (see 
§ 6.1). We then analyze three redshift bins 0.07 < 2 < 0.11, 
0.11 ^ z < 0.18 and 0.18 ^ z < 0.2, chosen to maximize 
the separation between high and low-z systems, while still 
retaining a significant number of clusters in each group. 
Fig. 12 presents the LFs derived for the three redshift 



Table 3. Results of the ML fit (within a projected radius of 
0.5 Mpc splitting our sample by redshift at z = 0.10 and z = 0.17 
for both the richest (6 ^ Richn.ML) and poorest (4 ^ Richn.ML < 
6) clusters. 
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bins. We find a clear indication that the dwarf to giant 
ratio increases with decreasing redshift, both for rich and 
poor systems. This confirms results (of varying significance) 
already reported in literature (Kodama ct al. 2004; Goto 
et al. 2005; Tanaka et al. 2005; Stott et al. 2007; Crawford 
et al. 2009) for cluster samples spanning much wider red- 
shift ranges. Similar trends have also been reported for field 
galaxies (Willmer et al. 2006; Xia et al. 2006; Ryan ct al. 
2007). We underline that the above result might hv affected 
by the decrease in completeness limit as a function of red- 
shift. We measure this effect by fitting data in the lowest 
redshift bin, using the completeness limit of the highest bin. 
We indeed observe a decrease in the steepness of the faint 
end (a = — 1.06l|5'()5 ^' ^ Richn.ML), but still signifi- 
cantly different from what measured in the highest redshift 
objects of our sample (see Table 3). It is also true that the 
cluster catalog completeness, as a function of redshift, must 
be taken into account in order to draw conclusions about the 
LF evolution with lookback time. The completeness of the 
NoSOCS catalog has been tested through extensive mock 
cluster simulations, as discussed in Gal et al. (2009): the cat- 
alogue is > 80% complete over the redshift range explored 
here for rich clusters, while for poor systems the complete- 
ness is a strong function of redshift, dropping below 50% 
at > 1.5 (see e.g. Figs. 4, 5 and 6 of Gal et al. 2003, for 
details). Finding the same redshift evolution in a in both 
rich and poor clusters thus strengthens our conclusions since 
the former are less subject to completeness bias due to the 
NoSOCS catalog, than the latter ones. Furthermore, any 
NoSOCS completeness effect would result in a steepening of 
the faint end at high redshift, since the sample there would 
be dominated by rich systems which have a steeper faint-end 
slope than rich ones (Figures 10 and 11). In this respect our 
result must be considered a conservative estimate of the LF 
dependence on lookback time. 

We further measure a trend for stronger suppression 
of faint galaxies (below M* -|- 2) with increasing redshift in 
poor systems, with respect to more massive ones, indicat- 
ing that the evolutionary stage of less massive galaxies de- 
pends more critically on the environment. A similar trend 
has been observed for faint red galaxies by Koyama et al. 
(2007), while discordant results are instead found by Craw- 
ford et al. (2009); Andreon (2008). 
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Figure 10. ML fit of the LFs extracted within 0.5 Mpc, 1.5 Mpc, and for 0.5 Mpc < R ^ 1.5 Mpc in bins of increasing richness. Relative 
confidence levels are also plotted. Results are for Richn.ML richness estimates. Data points show the composite LFs obtained using the 
GMA approach. 



Table 2. Results of the ML fit of the LFs extracted with different extraction radii in bins of increasing 
richness. Results are for Richn.ML richness estimates. 

Richn.ML < 4 4 ^ RichnML < 6 6 ^ RichnML < 7 7 s£ RichnML < 9 9 s£ RichnML 

R^ 0.5 Mpc 

a -U.01_goi -U-oa_o.o4 -^■'^^-i).o4 -l-U-i-o.os -l-l'J-o.o4 

M* -20.90l°;i2 -21.17l°;°5 -21.23l«:»^ -^^-^^it^o.ll -^^-^^to.it 

R^ 1.5 Mpc 

„ _nq+0-2 _1 ir;+0.08 _1 17+0 05 24+O.O6 _12^+°04 

a U.»_g 4 1.10_Q Qy ^-^'-0.05 ^-^^-O.OS ^-^"^-0.04 

M* -21.0+1-1 -21-47i°;l^ -21-39l°:i^ -21-46l°;i^ -21.56l°;l° 

0.5 Mpc s: R sS 1.5 Mpc 

_1 99+0-21 , Qn+0 06 _i 9f,+0.05 9Q+O.O4 _-, 90+" "4 

M* -21A+°i, -21.5l°;2 -21.37l°;i5 -21.34l0;il -21.54l0;i^ 

R< 0.5 R200 

_0 2+°-2 -0 92+° °'' -0 80+°-^° -0 95+° °^ -1 22+° °'^ 

M* -20.9l°j -21.28l°;J^ -21.29to,ll -'^^■^^to.H -21.48l°J^ 

R^ R200 

n «+0-3 1 Qo+0.05 -I nfi+0 07 , ns+0-05 i 00+0.04 

M* -21.21°^ -21.46l°:13 -21.29l°l? -21.25l°;° -21.30l°°« 

0.5 R200 s;r^ R200 



_1 0+0-3 _-, ,-,+0.10 _-, 99+0-08 _-, 15+0-10 _-, -,g+0-09 

« l-U_o.4 ^-^^-0-10 ^-^^-0-10 ^•^■^-0-08 ^-^"-O-OQ 

M* -20.91?* -21.21°:^ -21.2l0j -21.llgJ -21.3l°J 



The Luminosity Function of the NoSOCS Galaxy Cluster Sample 13 



0.5 R200 R200 




Magr 



Figure 11. ML fit of the LFs extracted within 0.5 R200, ^200 and for 0.5 R200 < R ^ R200 in bins of increasing richness. Relative 
confidence levels are also plotted. Results are for Richn.ML richness estimates. Data points show the composite LFs obtained using the 
GMA approach. 
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Figure 12. Left panel: ML fit of the LFs (extracted within 0.5 Mpc) of the richest structures (Richn.ML > 6), split into three redshift 
bins, together with relative contours, corresponding to the significance levels of 1 2 and 3 cr. Right panel: same as left panel, but for the 
poorest structures (Richn.ML < 6). 



7 CONCLUSIONS 



We presented the analysis of the Luminosity Function of 
galaxy clusters from the Northern Sky Optical Cluster Sur- 
vey, using r'-band data from the Sloan Digital Sky Survey. 
The sample, including 1451 galaxy groups and clusters, is 
large enough to allow us to investigate in detail both the 
intrinsic differences in the galaxy populations as a function 
of richness, cluster-centric distance and redshift, as well as 



the uncertainties introduced by different analysis techniques 
commonly used in the literature. 

Our global LF agrees with previous studies of galaxy 
clusters and does not show the presence of an "upturn" at 
faint magnitudes down to Mr ^ —16, presented by some 
earlier works as the proof of the presence of a very large pop- 
ulation of dwarf galaxies (e.g. Popesso et al. 2005; Gonzalez 
et al. 2006). 

We do observe a strong dependence of the LF shape 
(M* and faint-end slope) on both richness and extraction rar- 
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dius as expected from the morphology-density relation and 
most physical models of galaxy evolution in dense/massive 
systems. The dwarf to giant ratio increases with richness, 
indicating that more massive systems are more efficient in 
creating/retaining a population of dwarf satellites. Further- 
more, in the innermost regions R ^ 0.5 R200 we observe 
a sharp steepening of the faint-end and an M* brightening 
with richness. The same effect is observed both within fixed 
(0.5 Mpc) and physical (0.5 R2()()) apertures, suggesting that 
either the trend is due to a global effect, or to a local one 
but operating on even smaller scales, and thus giving rise to 
similar effects. Outside this radius, cluster galaxies appear 
to share similar LFs regardless of the mass of the object. The 
general trend for the LF to become shallower with decreas- 
ing cluster-center radii supports the hypothesis that dwarf 
galaxies are tidally disrupted near the cluster center, hence 
providing strong evidence that the relative mixture of giant 
and dwarf galaxies depends on the fraction of the virial ra- 
dius that is explored. This also explains why some studies 
using large aperture radii, even if physical, have missed these 
trends in the past. Our data further suggest that a signifi- 
cant growth of this dwarf population has occurred at rela- 
tively low rcdshift {z < 0.2) both in rich and poor systems. 
Both the richness and radial dependence of the faint-end 
slope are likely due to different mixtures of red/passive and 
blue/starforming galaxy populations, as observed by several 
authors (i.e. Zandivarez et al. 2006; Barkhouse et al. 2007), 
but the dwarf-to-giant ratio of red and blue galaxies has not 
been investigated in detail here and should be addressed in 
a forthcoming work. 

We note that an appropriate richness definition is re- 
quired if we want to extract information about the environ- 
mental effects on galaxy evolution based solely on optical 
data, since galaxy counts alone spanning a large magnitude 
range will introduce correlations between the LF slope and 
richness, which tend to dilute the observed differences. 

The results of LF studies such as the one presented 
here may also depend strongly on the input cluster catalog. 
Gal et al. (2009) compared the NoSOCS cluster catalog to 
the one derived from SDSS using the MaxBCG algorithm 
(Koestcr ct al. 2007). Due to the bright magnitude limit of 
DPOSS, NoSOCS is an essentially flux-limited sample with a 
richness-dependent completeness even at ^; ~ 0.2. As shown 
in fig. 3 of Gal et al. (2003), at highest richness the percent- 
age of NoSOCS recovered clusters is expected to be essen- 
tially redshift independent down to 2 ~ 0.2, while for the 
less rich groups analyzed in the present study {Ngai ~ 25; 
see Sec. 2.1) the percentage of recovery is expected to de- 
crease by a factor of ~ 2 between z ~ 0.07 and z ~ 0.2, 
with contamination rate being still smaller than ~ 5% (see 
fig. 8 of Gal et al. 2009). In contrast, the MaxBCG method 
relies on the E/SO ridge-line to detect clusters, and sam- 
ples such galaxies down to 0.4L* out to 2 = 0.4. Thus, 
the MaxBCG catalog, trimmed to 2: = 0.3 to reduce pho- 
tometric redshift uncertainties, provides something close to 
a volume-limited sample. Nevertheless, the requirement of 
a recognizable E/SO ridge- line may favor systems with un- 
evolved galaxy populations, resulting in different LF trends 
than those observed here. Indeed, (Gal ct al. 2009) find 
that at low richness both MaxBCG and NoSOCS may likely 
miss a significant firaction of groups. In fact, the MaxBCG 
catalog, restricted to systems with Ngais,MaxBCG > 10, 



misses about half of the poorest NoSOCS systems, de- 
spite the < 5% contamination of NoSOCS; on the other 
hand, NoSOCS also misses many of MaxBCG clusters with 
15 > Ngais, MaxBCG > 10, as cxpccted from the incomplete- 
ness of NoSOCS in this low redshift regime. This shows how 
the comparison of LFs of low-mass systems remains a chal- 
lenging issue, mainly because of the different selection func- 
tion and detection strategy utilized to construct different 
cluster catalogs. Thus, future studies of this population will 
require joining cluster catalogs generated using different al- 
gorithms to attain a fuller picture of the true underlying 
population(s) and the intrinsic LF variations. 

Finally we point out that LF studies based on large 
samples of different S/N clusters and groups are extremely 
sensitive to the technique used to both sum galaxies and to 
fit the galaxy distributions since the LF is far from universal. 
The uncertainties introduced by the different methods may 
explain in part the variety of faint-end slopes reported in the 
literature as well as, in some cases, the presence of a faint- 
end upturn. It is clear that the only way to prevent such 
uncertainties would be to use data probing the same abso- 
lute magnitude range for all clusters, such as what could be 
provided by the next generation of large surveys. Otherwise 
extreme care must be used in evaluating the effects that sta- 
tistical methods, in addition to the standard measurement 
errors, have on the final results. 
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